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Abstract 

A novel method for controlling nanowire magnetic properties and growth from filling time profile 
is presented. 

The wires are grown with an electrodeposition method (" Template synthesis" ) with a wide selection 
of pore diameters. We show that stray-fields presence in ferromagnetic nanowires are entirely 
dependent on the nanowire diameter. Besides a crossover effect in the reversal mechanism is 
observed with change in diameter. In this work, theory and experiment agree and confirm that 
according to the variety of hysteresis loop measured, about four ranges of values of pore diameter 
control the orientation of nanowire magnetization easy axis with respect to the geometrical axis. 
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I. INTRODUCTION 

There has been an increasing interest [lj in the fabrication and properties of nanostruc- 
tured magnetic materials not only from a fundamental but also from a technological point of 
view, owing to their potential applications in magnetic recording media [2J, biosensors and 
magnetic sensors based on the giant magneto- resistance effect 3|, |4| . The ultimate density 
of recording media depends on the size of the individual magnetic elements. The minia- 
turization of the latter can now be realized through the electrodeposition of nanodots or 
nanowires into self-assembled arrays, such as track-etched polymer membranes and anodic 
alumina filters. 

Structural and magnetic characterization of arrays of Nickel (Ni) nanowires produced by 
electrodeposition in polycarbonate membranes with diameters in the (15-100 nm) range and 
L = 6fim average length are made. The nanowires made within polycarbonate membranes, 
with various diameters (15 nm, 50 nm, 80 nm, 100 nm), are used to study the influence of the 
various parameters such as pore geometry and deposition process on the magnetic properties. 
Magnetization curves with various magnetic field orientations and nanowire diameters were 
determined at room temperature by magnetometry. According to the form of hysteresis 
loop measured, about four ranges of values of pore diameter control the orientation of easy 
axis of magnetization with respect to the wire axis. Reducing the diameter of the nanowires 
from 100 to 15 nm leads to increasing coercive fields from 347 to 590 Oe. The measured 
coercivity as a function of angle (y?) between the field and wire axis reveals that the coercive 
field decreases (increases) with angle, peaking at tp = 0°, for nanowire diameters smaller 
(larger) than 50nm. 

"Template synthesis" is an elegant chemical approach to the manufacturing of nanos- 
tructured materials, in particular for different kinds of nanowires [5j. In this paper, arrays 
of Ni nanowires have been obtained by filling a porous polycarbonate membrane, wich con- 
tains a large number of cylindical holes with a narrow size distribution. Characterization 
and understanding of the magnetic properties of nanowire arrays have been a challenge for 
years. Some problems remain unclear. For example, there are still open questions about the 
mechanisms responsible for the magnetization reversal. The intrinsic propeties of nanowire 
arrays are directly related to the properties of the nanoporous template such as the relative 
pore orientations in the assembly, pore size and its distribution, as well as interpore distance. 
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Two reversal modes have been suggested as being important: curling and coherent rota- 
tion. The critical diameter between coherent rotation and curling is d c = 2y^-^-, (d c is also 
called the coherent diameter), A = 1.5 x 10~ 6 erg/cm is the exchange stiffness constant 
Ms the saturation magnetization and q is a solution of a Bessel equation p| (see section 
3.1). 

For wire diameter d < d c , the Stoner-Wohlfarth coherent rotation model applies 

On the other hand for wire diameter d > d c , the reversal occurs by inhomogenous reversal 
(curling). Increasing further the diameter until d 3> d c , domains may form within the wire 
and the magnetization reversal, thus, may occur by domain wall motion 

As other extrinsic or hysteretic properties, coercivity (or coercive field) He, is strongly 
real-structure dependent. It is one of the most important properties of magnetic materials 
from the viewpoint of their technological utilization [h]] and the understanding of its intrinsic 
mechanism is a permanent challenge. As different magnetization reversal mechanisms would 
give a different angular dependence of the coercivity He, the measurement of Hc(ip) would 
provide helpful information about the rotation mechanisms. Here, <p is defined as the angle 
between the field direction and the wire axis (Fig.l). For anisotropic, perfectly oriented 
polycrystalline magnets based on noninteracting particles, the coercive field He should be 
equal to the anisotropy field Ha- It is known [ll|, however, that the coercive field measured 
in permanent magnet materials is typically one order of magnitude lower than the anisotropy 
field. It is basically accepted that the origin of this discrepancy lies in the distribution of local 
alterations of the magnetic properties due to a variety of microstructural features giving rise 
to magnetization reversal through a process wich includes a succession of mechanisms, as well 
as in the role played by interactions, either dipolar and/or exchange, during demagnetization. 

Aside from coherent rotation, important coercivity mechanisms are, curling, localized 
nucleation, and domain-wall pinning. With increasing size, the nucleation mechanism in 
perfect ellipsoids of revolution changes from coherent rotation to curling. However, both 
coherent rotation and curling greatly overestimate the coercivity of most magnetic materials. 
Nucleation refers to the onset of magnetization reversal and determines the coercivity in 
nearly defect-free magnets. 

Most works have concentrated on individual wires Q, few Q, flo1 | reported on 
detailed measurements of Hc(<p) for nanowire arrays. 

In this work, we set out to study the magnetic properties of various types of nanowires 
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possessing different diameters. Nanowires used in this work are fabricated in polycarbonate 
membranes, with various diameters (15 nm, 50 nm, 80 nm and 100 nm) in order to study 
the influence of the parameters such as pore shape and deposition process on the magnetic 
properties. 

We infer from our study that basically four ranges of diameter values are important to 
determine the orientation of the easy axis with respect to the nanowire axis. 



II. EXPERIMENTAL 



A. Sample preparation 



The nanowire samples are prepared by electrodeposition into the pores of nuclear-track- 
etched polycarbonate membranes PCTE (see Fig. 2) using an Autolab Potentiostat em- 
ploying a three-electrode configuration and an electrolyte contaning ion of the metals to be 



deposited [16] . 



During sputtering, a Copper layer was sputter-deposited on one side of the polycarbonate 
membrane and used as the working electrode to fabricate an array of Ni nanowires. A carbon 
plate and Ag/AgCl are used as the counter electrode and reference electrode, respectively. 
The work was carried out at a constant potential —IV (Ag/AgCl) at room temperature. 
The following aqueous electrolytes used was Nickel sulfate salt (NiS04.H 2 0)(44g/l) in a 
boric acid solution (H 3 B0 3 )(40g/l). The latter was used as a buffer. The study of the 
magnetic properties of the wire demanded that deposition should stop when there is no 
covering of the upper surface of the membrane by a polluting layer of the magnetic material. 
During the deposition process, the time dependence of electrical current was monitored and 
recorded. Hence, the current intensity-time curve recorded during the electrodeposition 
process revealed four different regions (see Fig. 3): 

- Region I: coverage of the pore walls by the ions. 



- Region II: filling of the pore interior by the growth. 



- Region III: fast variation triggered by growth at the upper extremity of the pore and 
subsequent outgrowth beyond the pore. This occurs when the pores are completely filled 
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with the material, and the electrodeposited material begins to form hemispherical caps over 
the nanowire ends (see Fig. 4). 

- Region IV is percolative growth outside the pores. 

Attention has been focused on the region labeled II because it is always the most com- 
plicated to reproduce. 

To understand the mechanism of growth materials in porous membrane, and to estimate 
the time of filling in the pores, a small comparative study between two electrodes was 
made. We tack two Copper electrodes in the same area, which was stuck on the first porous 
membrane (electrodel). According to the nickel growth curves on two electrodes (current 
versus time), we can clearly see the difference between the two curves because on the first 
membrane electrode (electrode 1), the deposit is made on small surfaces while on the another 
takes place over the entire surface. After some time, the two curves intersect in a single point, 
which shows that the pores are filled by materials and it is time for the emergence of a film 
on the surface of the membrane (Region IV, Fig. 3). This also indicates that deposition 
occurs on two equal surfaces (Fig. 5a). 

Figure 5b shows an image of a sample in which all the nanowires are misaligned per- 
pendicularly to the polycarbonate membrane plane. The majority of nickel nanowires tilt 
randomly a few degrees, with respect to the polycarbonate membrane plane normal. 

B. Magnetization measurements 

The magnetic properties of the Ni nanowire array are studied by VSM (vibrating sample 
magnetometer), measurements which are performed at room temperature. The hysteresis 
loops at different angles between the magnetic field and the axis of the wires were measured. 

Figure 6 shows four cases for the hysteresis loops with the external field H, perpendicular 
and parallel to the sample plane (i.e. parallel and perpendicular to the axis of the nanowires) 
for 4 samples with different diameters of (a) 15, (b) 50, (c) 80 and (d) 100 nm respectively. 

The hysteresis loops and the coercivities change smoothly with the change in the applied 
field direction. 

The saturation field, when the field is applied parallel to the axis of the nanowire, is 
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greater than that when the applied field is perpendiculary to the nanowire for the 15 nm 
diameter sample. On the other hand we see the opposite case for 80 and 100 nm diameter 
samples. For the 50 nm sample the parallel and perpendicular saturation fields are equal. 
(See Fig. 6). 

The coercivity dependence on the angles is shown in (Fig. 7a) demonstrating that co- 
ercivity decreases with the angle for nanowires that have 15 and 50 nm in diameter. In 
contrast, this is not the case of nanowires with 80 and 100 nm diameters. 

FIG. 7 (a-b) show the angular dependence of Hc{f) and SQ(<p) = j^- (squareness) for 
Ni nanowire arrays in PCTE membranes, respectively. 

However, Hc{f) shows a quite different behavior for, d=15 nm, 50 nm and d=80 nm, 
100 nm (Fig. 7a). For <i=15 nm and 50 nm, we observe a bell-type Hc((fi) variations, but 
for d =80 nm and 100 nm a bowl shape Hc((p) curve is observed. 

These results imply that there is a crossover transition in magnetization reversal mecha- 
nisms for d < 50 nm. 

As shown in (Fig. 7b), for the two samples, 80 and 100 nm diameters, the squareness 
SQ(ip) of Ni nanowire array decreases with increase of angle, peaking at ip = 0°, which is 
the easy axis of the arrays as can be seen from the magnetization curves (Fig. 6). 

On the other hand, for the sample of 15 nm diameter the easy axis is perpendicular to 
the wire axis. For the 50 nm diameter sample, easy magnetization axis is off to an angle of 
30° of the wire axis. 



III. THEORETICAL 



In this work we can infer that single domain in the nanowires, and magnetization reversal 
modes can be modeled by homogeneous rotation when there is a critical size below which 
a particle remains in a single-domain state during switching, or inhomogeneous space de- 
pendent reversal 8(r) when the particle size is larger then the critical size, but still in the 
single-domain regime. We assume that all our samples are in the single-domain regime 
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A. Magnetization in single-domain particles 



In principle, the magnetization configuration in a magnetic nanowire can be determined 



from the Brown equations by minimizing the total free energy [171 ] . In bulk ferromagnetic 
materials, the energy of the system can be minimized by forming multiple magnetic domains 
within which the magnetic moments are aligned. However, there is a critical size below which 
a particle remains in a single-domain state during switching. 

Approximating our single domain by an ellipsoid of revolution, the critical radius R s d 
(minor axis, a = b) for a single-domain particle can be expressed and calculated by 
comparing the exchange energy averaged over the ellipsoid volume to the magnetostatic 
energy; [Appendix A]. One finds: 



T3 AR sd 

Where Nq is the demagnetization factor along the c axis (Fig. 1) and ai is the nearest- 
neighbor spacing. The critical radius for a single-domain particle is dependent on the ma- 
terial parameters A, a\, Ms- 

If we consider the infinite cylindre case (R s d —* oo,Nc — > 0), in the single domain case 
there is a critical size between the two magnetization reversal processes. When the magnetic 
ea 8y a* is afigned with the apphed field, the critical radius r c = f for the transition between 
coherent and curling effect is defined by [7j: 

l2A q 

where q = 1.8412 for a cylinder is the first maximum of the Bessel equation -^-(Ji(q)) = 
0, where J\(q) is the ordinary. Bessel's function of the first kind and N a = 2ir is the 
demagnetizing factor along the minor axis of the infinitely long cylinder. Hence: 

^4 q 



7T M S 

Coherent and curling model rotation are illustrated in Fig. 8 

For the theoretical explanation of these two magnetization modes, we consider the elon- 
gated single-domain particle with uniaxial anisotropy shown in Fig.l. 
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B. Coherent rotation 



In the Stoner-Wohlfarth (SW) coherent rotation model 18J, the total free energy consists 
of the crystalline and shape anisotropy and Zeeman energy due to the external magnetic 
field. 



E = Ksm 2 9-M s Hcos{9-cp), (1) 

Where (9) is the angle between the Ms and the easy axis, K is the effective anisotropy, 
(<p) is the angle between the external applied field and the easy direction; {9 — (p) is the 
angle between Ms and H. 

M 2 

K=^(N a -N c ) + K 1 

Here N a and Nq are the demagnetizing factor of the ellipsoid along the minor and the 
long axis respectively, and K\ is the magnetocrystalline anisotropy contant, K\ is small by 

JUT 2 

comparing with -^-{N a — N c ). 

After using f§ = and = we obtain: 



He 



2 2 3 

[siny)3 + cosyjsji 



where Hs is the switching field and Hk is the effective anisotropy field : Hk ~ j/v 



M S - 



(The detailed calculation is in appendix B). 



On the other hand, the coherent rotation mode 17[ gives the highest and the lowest 
coercivity field Hc{<p) value of H parallel and perpendicular to the easy axis, respectively, 
for the 15 and 50 nm of diameter. 

In the Stoner-Wohlfarth model the switching field ( Hs) , does not represents the coercivity, 



Hc(<p), in all cases. However, from the discussion in 18] ; the coercivity can be written as 



, Hex 10 < <p < tt/4, 
H c = { fori 

H C 2 vr/4 < if < tt/2, 



with 



(2) 



Hci = \H& 



2K 



1 



M : 



s [sin cp 3 + cos tp s J 2 



H, 



C2 



Hs(f = t) 



Hs(\ < „ < £ 



In this mode, two regimes can be identifed, when the applied field is close to the magnetic 
easy axis, 0° < tp < 45°, the hysteresis loop is relatively square; hence, the change in sign of 
the magnetization, corresponding to the coercivity, occurs at the switching field. 
In this regime the coercivity field Hex is equal to the switching field H$- 
For 45° < <p < 90° , the applied field is oriented closer to the magnetic hard axis, and the 
hysteresis loop is sheared such that switching occurs after the magnetization changes sign. 
In this case, He ^ Hs and the coercive field is determined from equation (Hc2)- 



C. Curling 



The magnetization curling mode was defined by Frei et al 17| and after it has been used 



for different strnctnres to investigate the magnetic switching of films 19], spherical particles, 
prolate ellipsoids [17[ and cylinders [20]. For particle size larger than the critical size but 
still in the single-domain regime, magnetization reversal occurs by curling {itI \\\ . In the 
curling model, magnetization switching is an abrupt process, and the switching field is very 
close to the nucleation field; hence, H c = H$ for all angles. 
The curling nucleation field is given by |21[ ] . 



h cos(^ -9) = ^-(N a sin 2 9 + N c cos 2 9) - 

Z7T D 



(3) 



where S is the reduced radius defined as R/Ro, R is the radius of the cylinder and R0 is 
the exchange length defined by Rq = and Rq = 25.2 nm for Ni. Here h = 2 Jif s 1S the 
20 | and the parameter k = 1.08, which is q 2 /n in the notation of Ref. [3], is a 



reduced field 



monotonically decreasing function of the aspect ratio m = - of the ellipsoid (Fig. 1). 
To complete the calculation it is still necessary to eliminate the angle 9 from (3). 
By using the fact that it is the same angle as in: 



hsm(<f -9) = ^-(N a - N c ) sin2f9 
2n 



(4) 



it will become after the differentiating E with respect to 9. 
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dE 

— = 2K sin 9 cos 9 + M S H sm(9 -<p)=0. 

o9 

One way of doing it to consider Eqs (3) and (4) as two linear equations in h sin <p and h cos p 
and to solve them as such obtaining: 

1 k 
hsinp = (^N a - — ) sin# 

1 k 

hcosp = (—N c - —)cos9 (5) 
Adding the squares of these equations by using: 

cos 2 9 + sin 2 9 = 1 

leads to 

hs 



( N c k_\(Na k_\ 



V(|?-^) 2 8inV + (^-^) 2 co8V 

Wllere ^ = ^fe- 

For the infinite cylinder, Nc = 0, N a = 2n, and = 1.08 we obtain 

1.08(1 - 1.08g- 2 ) 

S ~ C ~ ,S 2 [1.1664^ 4 + cosV(l -2.16^ 2 )]5' 

The coercivities of infinite cylinders with shape anisotropy depends only on the value of 
S and the angle tp between the easy axis of magnetization and the direction of measurment 
as shown in (Fig. 9-10) and plotted in reduced units h c (ip) (0 < h c < 1) for various values 
of S. H c versus cp curves are illustrated in Fig. 11. When we compare the curves in Fig. 7a 
and Fig. 11, the coherent magnetization reversal is for the 15 nm and 50 nm of diameters 
and the curling mode is for the 80 nm and 100 nm of diameters. 

In order to observe small variation in the coercivity, very well geometrically characterized 
samples need to be measured with well controlled inter-wire interaction. In our samples the 
center to center distance is fixed, and the strength of interactions is different from sample 
to sample, making direct comparison difficult. The fitted values for the coercivity are larger 
than the experimental data and Hc(<p) values are actually much smaller than 2n Ms = 3050 
Oe, which is expected for an individual wire. 
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We ascribe such differences between calculations and experimental results to the interac- 
tion of each wire with the stray field produced by the array. This field originating from the 



effective ferromagnetic coupling between neighboring wires reduces the coercive field [221.123] . 

Fig. 6 a, shows there is a flat behaviour in the Hc(<p) curve at low angles for 80 and 100 
nm of diameters, which may be caused by the variation of wire orientation with respect to 
the normal to the membrane plane. 

These results suggest that the interwire coupling cannot be neglected even for PCTE 
membranes. 



IV. DISCUSSION 

In the case of a single domain wire, we have to consider three different field acting on the 
wire [24]: 

Heff — H S h + H A + H D 

where Hsh is the demagnetization field (shape), Ha is the effective anisotropy field and 
Hd is the average dipolar field between the wires. 

The demagnetization field (shape) H$h of the individual wire if magnetized parallel to 
the pore axis is of the order of 2nMs = 3050 Oe with Ms equal to 485 emu/cm 3 . 

The effective anisotropy field Ha, given by Ha = .~f = 123 Oe with K\ = —4.5 x 10 4 



erg/cm 3 is the cubic magnetocrystalline anisotropy [25| of Ni. 

The in plane average dipolar field between the wires, which tends to induce a magnetic 
easy axis perpendicular to the wire axis. 

For a qualitative analysis, we consider a two-dimensional infinite array of magnetic 
dipoles. The array is formed by a square dipole mesh of size D. When we assume all 
the moments are aligned along the wires, the average dipolar field, which acts on each wire, 
is parallel to the axis of the wire and may be expressed as 



2j: 



H Do = [A.2M s (nd 2 /A)L]/D 3 (6) 

On the other hand, for the same array with magnetic moments aligned perpendicularly 
to the axes of the wires, the average dipolar field which acts on every wire may be expressed 



as 



26. 



23| : 
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H Dgo = [~2.1M s (nd 2 /4)L]/D 3 (7) 

where d is the diameter of the wire, L is the length, and D is the inter-wire distance 
[Appendix C]. 

While the total effective fields following the two directions parallel (0°) and perpendicular 
(90°) to the wire axis may be written: 

' (0°) -> H Do + H A 

H e ff (total) = < 

(90°) -> H D90 + H Sh 

The calculation for the various angles, show that Hc(f) decreases as the angle increases, 
with a bell-type variation. Therefore, the magnetization reversal was determined mainly 
by the curling mode, as shown in (Fig. 6a) for 80 nm and 100 nm diameters. For Ni 
nanowires, the curling is predicted for d > 50 nm, this value agrees with prior work on 



ordered Ni nanowire arrays in alumina 23| where the critical diameter was 40 nm. The 
coherent rotation will occur for 15 and 50 nm wires, which is consistent with results in (Fig. 
6a). In the same way, the squareness SQ((p) of hysteresis loops is measured for different 
orientations of wires as a function of the applied magnetic field. It can be seen that for large 
wire diameters the maximum SQ(ip) is found when the magnetic field is applied parallel 
to the axis of the wire and the minimum occurs when the applied field is perpendicular to 
the axis of the wire (see Fig. 6.b). This concurs with the abovementioned idea that for the 
larger wire diameters, the easy axis of magnetization is parallel to the axis of the wires. On 
the other hand, for the 50 nm of diameters of wires (see Fig. 6.b), the behavior of SQ(ip) 
indicates that we are in the crossover region where the axes begin to rotate from parallel to 
perpendicular to the wire axis. But in the wires of 15 nm of diameter, the behavior of SQ(<p) 
indicates that we are in the perpendicular to the axes of the wires becomes the dominant 
easy axis of effective anisotropy. 



V. CONCLUSION 



In conclusion, by means of theoretical studies and experimental measurements, we have 
investigated the reversal processes in ferromagnetic nanowires. Our systematic studies of 
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the effect of the nanowire size diameter show that the magnetization reversal mechanism is 
strongly influenced by the diameter value of the nanowires. The various types of Hc((p) 
curves observed in (Fig. 7) suggest that no simple magnetization reversal mode could 
account for the complex coercivity mechanism. Two reversal modes are considered as the 
most important: coherent rotation and curling. Good agreement between the measured 
magnetic properties of Ni nanowires and the theoretical calculations is obtained. However, 
further experimental work remains to be done in order to observe this transition. 
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VI. APPENDIX A 

Critical radius R sd for single domain 

We consider in Fig. 1, a — b — R as the short axis, and c as the long axis of the ellipsoid 



dr r l 



Approximating: 0(r) 



| The average over ellipsoid: 




with equation: ^2+^- = l^>Z = ±cJ 1 





where a\ is the nearest-neighbor spacing 




13 



- 3A r R dr l 
Ex = —j a —c\\l- 



2 

R 2 



changing: = cos a; =^> ^ = 



E 



R 2 c i a V 2 r 
sin xdx 

f^ccos( a y 2 R) Rsin 2 x 



6A r a 
R 2 Jo 



Rcosx 



dx 



3A f^ccos( a y 2 ji) s i n 2 x 



_ 6A [ 
~~ R 2 Jo 



dx 



cos x 



From the integral: / ^fdx = —[sin a; + hitg(f - §)] 
Taking: u = arccos( a i/2i?) =>• a V2.R = COSM 

cos sin 77 



using: _B = \ntg(j — |) = In 
witht = t,f,5 = l T f=(l5# = ^ 



5 



COS M 



1 sin m 1 - V' 1 - (|r) z 



cosm 



ni 
2_R 



R - 3A 



1 + ln^i- 
4i? 



3A 
^2 



ln^-1 



Comparing this energy to the magnetostatic energy: 
We get: E D = \M 2 N C 



R 



sd 



r[ln(— ) - 1] 



N c M 2 s l ai 



VII. APPENDIX B 



Stoner-Wohlfarth (SW) coherent rotation model 
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The total free energy, (see Fig. 1) 

E = K sin 2 - M S H cos{9 - <p) 

After differentiating E with respect to 9, we obtain 

dE 

— = 2K sin 9 cos 9 + M S H sin(0 - <p) , 
o9 

We find that for a minimum in E, ^ = 0, i. e; 



2Ksm 9 cos 9 + HM S sin(# - <p) = 



After using the trigonometry equation: 

sin(6> — V 7 ) — sm $ cos V 9 ~~ cos $ s i n V 9 We find: 



r H Ms cos 09 i^M*? sin , 

2K+ A - y ^r~H sin#cos# = 

cos 9 sin 9 



H z = H cos ip 

Substituting: 

Hx = H sin (p 



If we exlude the extrema at 9 = 0, n and = |, , we can assume that above terms (in 



square brackets) is zero. 



// : .U S H X M S 

2K + —5—^ = 

cos 6 1 sm 9 



Defining: 2 2 * H « , H K = jf 

Ux 2K H K 

The second derivation gives: 



d 2 E h x sin 9 h z cos 9 
~W = cos 2 9 + sin 2 # 



After using fjf = and C-f = 0, we obtain: 



cos y sm v / 

/,; h-AB.2) 



Eliminating from these equation we obtain: 
By (B.2), ^ = 

Trough ( B .l),-gf/i,--J f » = l 
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Factorizing: ^p^[sin 2 9 + cos 2 9]h z — 1 
h z = — cos 3 9 



We get < 



h x = sin 3 9 



which gives us h x + /iz — 1, (Astroid equation) 

Replacing h x and h z by these values we find the switching field Hs 



2 2 3 

[sin<^3 + cos (ps] 2 



VIII. APPENDIX C 



Dipolar interaction field 



If we have a single dipole P surrounded by an array of dipoles p (Fig. CI), the interaction 

is: 

E _yP .P t 3(P .r)(P.r) 

i 

This is equivalent to a field H D acting on the dipole such that: E = —P .H D 

S ff ^ p o-Pi 3(P .f)(P.f) 
-P M D - ^ — 

i 

Calculating H Dz 

Assuming the dipoles Pj sit on a 2D square lattice of parameter D (the average nanowire 
separation) 



Pn 



hence P 



and P 



1, is directed along the Z axis, 


=> -P .H D = -H Dz 
1 

iD 

P , and f 



v 
Pz 



jD P .f=0,P t .P = P z , 


we have r 3 = D 3 (i 2 + j 2 )^ 
Then, H Dz = - E ^ 

Pz = M5T/(AU surrounded dipoles saturated along Z): 
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Dz ~ — W ^ 37 

id (i 2 + f) /2 



+00 

E — ^r- = 4 E 



,2 + J 2) 3 / 2 . _ (,2 + j2 ) 3 / 2 



i, j = —00 v ' % — 1, 00 



Numerically: 



4 X E !-3- ^ 4 x (1.08) 

l.oo < i2+ «'2 j = l lfl0 P+i^ 



J = l,oo j = 1,100 



-4.2M 5 y 



Calculating H Dx 

Taking, 

1 

Pn = , P ~- 




P Y , and r 
Pz 



%D 
JD , 




=S> E = -P .H D = -H Dx ,P .P = P x 



)E P x 3(iD)(P x iD + P Y jD) 



We get E = ^fi^ 



=>- = - 2^ 5 

ij=-oo £3^2 + J 2) /2 



If all the surrounding dipoles are saturated along X: 



P x = M s V,P Y = 
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wp „ P t ff n _M s V v oo (2i 2 -j 2 ) ^tt _ AM S V y^oo (2i 2 -f) 

we get n Dx - -^g- 2^ ij= „ 00 5— =^-«d x - -53- Z^ J= i 5— 

Numerically: 

± (2 ' 2 ~ i . ) = E ^0.526 

ij=l (j2 + j2 ) / 2 ij=l ( f 2 + j2 ) / 2 



^ X = 4(0.526)-A_ = 2.1-^1 



=>■ % - ~2 H d z 
This shows the dipolar field along x is half the field along Z. 
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Figures 
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FIG. 1: Single domain particle with uniaxial anisotropy in presence of applied field H, showing 
angles, 9 and ip that M and H make with wire axis along Z, (a = b < c). 
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FIG. 2: Polycarbonate PCTE membrane with 80 nm pore diameter. 



22 




Time (s) 

FIG. 3: Typical time variation of cathode current during electrodeposition of Ni nanowires arrays 
or we see the four regions, coverage of the pore walls by the ions (region I), filling of the pore 
interior by the growth (region II), fast variation trigged by growth at the upper extremity of the 
pore and subsequent outgrowth beyond the pore. This occurs when the pores are completely filled 
with the material, and the electrodeposited material begins to form hemispherical caps over the 
nanowire ends (region III), percolative growth outside the pores (region IV). 
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FIG. 4: SEM observation of the top surface where the Ni nanowires emerge from the polycarbonate 
PCTE membranepores in region III. 
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FIG. 5: Time variation of cathodic current during electrodeposition of Ni nanowires arrays when 
we can determine the exact time of filling of the pores (a) the two intensity curvus cut at point 
P where the vaalus of the top surfaces are the same. SEM image of a sample in which all the 
nanowires are misaligned perpendicularly to the polycarbonate membrane plane (b). 
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FIG. 6: a- The hysteresis loops when the external field is perpendicular and parallel to the axis 
wires for 15 nm and 6/im length. 

b- The hysteresis loops when the external field is perpendicular and parallel to the axis wires for 
50 nm and 6//m length. 

c- The hysteresis loops when the external field is perpendicular and parallel to the axis wires for 
80 nm and 6//m length. 

d- The hysteresis loops when the external field is perpendicular and parallel to the axis wires for 
100 nm and 6/j.m length. 
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FIG. 7: Angular dependence of the coercivity Hc{^>) (a) and the squareness SQ (b) of 
in PCTE membranes with different diameters: 15, 50, 80 and 100 nm and 6/im length. 
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FIG. 8: Illustration of magnetization reversal in a single-domain cylinder: (a) coherent rotation 
and (b) curling case. 
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FIG. 9: Angular dependence of the reduced coercivity, he = 2 ^M S ' * n Stoner Wolhfarth case, 
where <p is the angle between the external applied field and the wire axis, cZis the wire diameter, 
m = ^is the aspect ratio and S = j^is the reduced ratio. 



31 




FIG. 10: Angular dependence of the reduced coercivity, he = 2nM s ' * n ^ e Curling case, where <p 
is the angle between the external applied field and the wire axis, dis the wire diameter, m = -is 
the aspect ratio and S = -S^is the reduced ratio. 
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FIG. 11: Reduced coercivity h c for aligned array of infinitely long single-domain cylindrical 
nanowires as function of angle ip between easy axis (i.e. cylindrical axis) and direction of mea- 
surement, where ip is the angle between the external applied field and the wire axis, d\s the wire 
diameter, m = §is the aspect ratio and S = -S^is the reduced ratio. 
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FIG. 12: Dipole Po surrounded by a number of Pj dipoles sitting on a 2D square lattice with D 
the average dipole separation. 
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